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Abstract 

In this paper we show that uniformly global asymptotic stability for a family of 
ordinary differential equations is equivalent to uniformly global exponential stability 
under a suitable nonlinear change of variables. The same is shown for input-to-state 
stability and input-to-state exponential stability, and for input-to-state exponential 
stability and a nonlinear estimate. 
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1 Introduction 



Lyapunov's notion of (global) asymptotic stability of an equilibrium is a key concept in the 
qualitative theory of differential equations and nonlinear control. In general, a far stronger 
property is that of exponential stability, which requires decay estimates of the type || < 
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ce~ xt \\x(0) ||." (See for instance [16] for detailed discussions of the comparative roles of asymp- 
totic and exponential stability in control theory) In this paper, we show that, for differential 
equations evolving in finite-dimensional Euclidean spaces M, n (at least in spaces of dimensions 
7^ 4, 5) the two notions are one and the same under coordinate changes. 

Of course, one must define "coordinate change" with care, since under diffeomorphisms 
the character of the linearization at the equilibrium (which we take to be the origin) is 
invariant. However, if, in the spirit of both structural stability and the classical Hartman- 
Grobman Theorem (which, cf. [23], gives in essence a local version of our result in the 
special hyperbolic case), we relax the requirement that the change of variables be smooth at 
the origin, then all obstructions disappear. Thus, we ask that transformations be infinitely 
different iable except possibly at the origin, where they are just continuously differentiable. 
Their respective inverses are continuous globally, and infinitely differentiable away from the 
origin. 

Closely related to our work is the fact that all asymptotically stable linear systems are 
equivalent (in the sense just discussed) to x — —x; see e.g. [1] . The basic idea of the proof in [1] 
is based upon projections on the level sets of Lyapunov functions, which in the linear case of 
course be taken to be quadratic (and hence have ellipsoids as level sets). It is natural to use 
these ideas also in the general nonlinear case, and Wilson's paper [36], often cited in control 
theory, remarked that level sets of Lyapunov functions are always homotopically equivalent 
to spheres. Indeed, it is possible to obtain, in great generality, a change of coordinates 
rendering the system in normal form x = —x (and hence exponentially stable), and several 
partial versions of this fact have appeared in the literature, especially in the context of 
generalized notions of homogeneity for nonlinear systems; see for instance [6,25,15,27,24]. 

It is perhaps surprising that, at least for unperturbed systems, the full result seems not to 
have been observed before, as the proof is a fairly easy application of results from differential 
topology. (Those results are nontrivial, and are related to the generalized Poincare conjecture 
and cobordism theory; in fact, the reason that we only make an assertion for ^ 4, 5 is closely 
related to the fact that the original Poincare conjecture is still open.) 

Note, however, that it has been common practice in the papers treating the nonlinear case 
to use the flow generated by the original system to define an equivalence transformation, 
thereby reducing the regularity of the transformation to that of the system. Here we use 
the flow generated by the (normalized) Lyapunov function itself, which yields more regular 
transformations. In addition, and most importantly, our poof also allows for the treatment 
of perturbed systems (for which the reduction to x — —x makes no sense). 

Lyapunov's notion is the appropriate generalization of exponential stability to nonlinear 
differential equations. For systems with inputs, the notion of input to state stability (ISS) 
introduced in [29] and developed further in [5,9,13,14,17,18,26,28,32,33] and other references, 
has been proposed as a nonlinear generalization of the requirement of finite C 2 gain or, as 
often also termed because of the spectral characterizations valid for linear systems, "finite 
nonlinear H°° gain" (for which see e.g. [2,11,12,34]). We also show in this paper that under 
coordinate changes (now in both state and input variables), the two properties (ISS and 
finite H°° gain) coincide (again, assuming dimension ^ 4,5). 
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We do not wish to speculate about the implications of the material presented here. Obviously, 
there are no "practical" consequences, since finding a transformation into an exponentially 
stable system is no easier than establishing stability (via a Lyapunov function). Perhaps 
these remarks will be of some use in the further theoretical development of ISS and other 
stability questions. In any case, they serve to further justify the naturality of Lyapunov's 
ideas and of concepts derived from his work. 



2 Setup 

We consider the family of differential equations 



where / : R" x D -> R™ is continuous and for x 7^ locally Lipschitz continuous in x, where 
the local Lipschitz constants can be chosen uniformly in d G D C R m . Let T> denote the set 
of measurable, locally essentially bounded functions from R to D. For any xq G R n and any 
d(-) G V, there exists at least one maximal solution of (1) for t > 0, with x(0) = x . By abuse 
of notation, we denote any such solution, even if not unique, as (f>(t,x ,d(-)), t G I(x,d(-)), 
where I(x,d(-)) is its existence interval. Throughout the paper, || ■ || denotes the usual 
Euclidean norm, and "smooth" means C°°. For a differentiable function V : R n — > R the 
expression Lf d V(x) denotes the directional derivative DV(x)f(x, d). 

The general framework afforded by the model (1) allows us to treat simultaneously classical 
differential equations (the case when D = {0}) and more generally robust stability of dif- 
ferential equations subject to perturbations (when functions in V are seen as disturbances 
which do not change the equilibrium, as in parameter uncertainty), as well as systems with 
inputs in which elements of V are seen as exogenous tracking or regulation signals, or as 
actuator errors (in which case, the continuity properties of (x, d) \— > 0(-, x, d) are of interest). 
In light of these applications, we now describe the appropriate stability concepts. 

For the first, assume that D is compact and that /(0, d) — for all d G D. Then we say that 
the zero state is uniformly globally asymptotically stable (UGAS) if there exists a class ICC 
function (5 such that, for each d(-) G V, every maximal solution is defined for all t > and 



for all t > 0. As usual, we call a function a : [0, 00) — > [0, 00) of class /C, if it satisfies 
cn(0) = and is continuous and strictly increasing (and class /Coo if it is unbounded), and we 
call a continuous function (3 : [0, oo) 2 — > [0, 00) of class ICC, if it is decreasing to zero in the 
second and of class /C in the first argument. (It is an easy exercise, cf. e.g. [20], to verify that 
this definition is equivalent to the requirements of uniform stability and uniform attraction 
stated in u e — 5" terms.) Note that while our general assumptions on the right hand side / 
do not guarantee uniqueness of solutions through zero, the added assumption of asymptotic 
stability implies that 4>(t, 0, d) = is the unique solution with initial condition x — 0, for all 



x(t) = f(x(t),d(t)) 



(1) 



U(t,x,d(-))\\ <P(\\x\\,t) 



(2) 
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d G V. As a consequence, since away from zero we have a local Lipschitz condition, solutions 
are unique for each given initial state and d G V. 

If the origin is no common fixed point for all values d G D then (2) is impossible. In this 
case, however, still a useful notion of stability is possible. We call the system (1) (globally) 
input-to-state stable (ISS), if there exists a class /C£ function (3 and a class /Coo function a 
such that all solutions of (1) satisfy 

U(t,x,d(-))\\ < /Jfllxll.f) + «( sup || d ( T )||) (3) 

0<T<t 



for all d(-) G T> and all t > 0. Formulation (3) is the most frequently used characterization 
of the ISS property. Note that with (5 — 2(5 and a = 2a inequality (3) immediately implies 

\\(j>(t,x,d(-))\\ < max j/3(||:r||, i), a( sup \\d(r)\\)\ , 



hence this "max" formulation can be used as an equivalent characterization. 

Two apparently stronger formulations of these properties are obtained if we replace t) 
by ce _A *||x||, more precisely we call the zero position of (1) uniformly globally exponentially 
stable (UGES), if there exist constants c > 1, A > such that 

U{t,x,d{-))\\<ce- xt \\x\\ (4) 



holds for all d(-) G V and all t > 0, and we call the system input-to-state exponentially stable 
(ISES), if there exist a class /Coo function a and constants c > 1, A > such that 

\\(f)(t,x,d(-))\\ < max Jce" A< ||x||, a( sup ||d(r)||) > (5) 

[ 0<T<t J 



for all d(-) G T> and all t > 0. (As usual, these definitions use appropriate constants c, A > 0. 
In this paper, however, we will see that we can always work with "normalized" versions 
choosing c — 1, X — 1. For the (ISES) property we use the "max" formulation because it 
allows a further implication as stated in Theorem 5, below. Observe that (5) implies (3) with 

(5(\\x\\,t) = ce- A '||x||.) 

Extending the concepts in [1, p. 207] to our nonlinear setting, we will call a homeomorphism 

T : W 1 -> W 1 



a change of variables if T(0) = 0, T is C 1 on R™, and T is diffeomorphism on R n \ {0} (i.e., 
the restrictions of T and of T _1 to R n \ {0} are both smooth). Given a change of variables 
T and a system (1), we may consider the transformed system 

y(t) = f(y(t),d(t)), (6) 
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where, by definition, 

f(y,d) = DT(T- 1 (y))f(T' 1 (y),d). 

In other words, system (6) is obtained from the original system by means of the change of 
variables y = T(x). Observe that the new system again satisfies the general requirements: 
f(y, d) is continuous, and it is locally Lipschitz on x for x ^ 0, uniformly on d. 

It is our aim to show that for dimensions n ^ 4,5 the following assertions are true. Given a 
system of the form (1) satisfying (2) or (3), respectively, there exists a transformed system 
that satisfies (4) or (5), respectively. In this sense, global asymptotic stability is equivalent to 
global exponential stability under nonlinear changes of coordinates. Furthermore, one may 
obtain transformed systems where the constants defining the exponential stability property 
can be chosen to be the special values c = A = 1. 

Furthermore we show that if system (1) is ISES (5) with c = A = 1 then there exists a 
homeomorphism R : R m — > R m on the input space with R(0) = that is a diffeomorphism 
on R m \ {0} such that the transformed system with v = R(d) 

y(t)=f(x(t),v(t)), f(x,v) = f(x,R~\v)) (7) 

satisfies the following "L 2 to L 2 " nonlinear estimate: 

t t 
J \\<t>{s,xM-))fds<\\x\\ 2 + J \\v(s)\\ 2 ds. (8) 



Since (8) in turn implies ISS (by [31, Theorem 1]), we obtain equivalence between ISS and 
the nonlinear estimate (8) up to nonlinear changes of coordinates. 



3 Construction of the coordinate transformation 

The main tool for our construction of T is the use of an appropriate Lyapunov function V . 
In fact, we can obtain T for a whole class of functions as stated in the following proposition. 
Recall that a function V : R™ — > R is called positive definite if V(0) = and V(x) > for 
all x 7^ 0, and proper if the set {x | V(x) < (3} is bounded for each (3 > 0. 

The next result says in particular that any such function may look like ||x|| 2 under a co- 
ordinate change. This implies in particular that the level sets under coordinate change are 
spheres. It may therefore not come as a surprise that a basic ingredient of the proof is related 
to the question of whether level sets of Lyapunov functions in R™ are diffeomorphic to the 
sphere S n ~ 1 . This question is solved except for the two special cases of dimensions n = 4 
and n — 5, though in the case n = 5 it is at least known that the statement is true if only 
homeomorphisms are required. (For the case n = 4 this question is equivalent to the Poincare 
conjecture; see [36].) 
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Proposition 1 Let n 7^ 4 and let V : R n — > R &e a proper, positive definite C 1 function. 
Assume furthermore that V is smooth on R™ \ {0} with nonvanishing gradient. Then for each 
class /Cqo function 7 which is smooth on (0, 00) there exists a homeomorphism T : R n — > R n 
mt/i T(0) = such that 

V{y):=V{T-\y))= 1 {\\y\\). 
In particular this holds for 7(||y||) = \\y\\ 2 ■ 

Ifn 7^ 4, 5 then T can be chosen to be a diffeomorphism on R n \{0}. Furthermore, in this case 
there exists a class /C^ function 7 which is smooth on (0, 00) and satisfies 7(s)/7'(s) > s 
such that T is C 1 with DT(0) = 0. 

PROOF. For the function V the right hand side of the normed gradient flow 

. _ VV(x)' 

x ~ liw(*)ll 2 

is well defined and smooth for x 7^ 0. Denote the solutions by ip(t,x). Then V(ip(t, x)) = 
V(x)+t, and thus since V is proper and W(x) 7^ for x 7^ for a given initial value 
•0 is well defined for all t £ (— V(x), 00), thus also smooth (see e.g. [10, Corollary 4.1]). 

Fix c> 0. We define a map vr : R" \ {0} -> V~\c) by 

7r(a;) = ip(c — V(x), x) . 

Obviously ir is smooth, and since the gradient flow crosses each level set V~ 1 (a),a > 
exactly once it induces a diffeomorphism between each two level sets of V, which are C°° 
manifolds due to the fact that V is smooth away from the origin with nonvanishing gradient. 

Now observe that the properties of V imply that V~ 1 (c) is a homotopy sphere (cf. also 
[36, Discussion after Theorem 1.1]), which implies that V~ l (c) is diffeomorphic to § n_1 for 
n = 1, 2, 3 (see e.g. [22, Appendix] for n = 2, [7, Theorem 3.20] for n = 3; n = 1 is trivial). For 
n > 6 we can use the fact that the sublevel set {x e R n | V(x) < c} is a compact, connected 
smooth manifold with a simply connected boundary, which by [21, §9, Proposition A] implies 
that the sublevel set is diffeomorphic to the unit disc D n , hence V -1 (c) is diffeomorphic to 
§n-i_ Thus for all dimensions n 7^ 4, 5 we may choose a diffeomorphism S : V _1 (c) — > § n_1 . 
By [8] we may choose S to be at least a homeomorphism in the case n = 5. 

Let Q := S o n. The coordinate transformation T is now given by T(0) = and 

T(x)= 1 -\V(x))Q(x), x^0. 
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An easy computation verifies that T 1 (0) = and 



T- 1 (y)=il>h(y)-c,S- 1 ( A , 2/7^0, 



hence T is a diffeomorphism on IR n \ {0} (resp. a homeomorphism if n — 5). Since ^(0) = 0, 
and ■?/>(£, — > as t \ — c, both T and T _1 are homeomorphisms. 

Finally, we have that 



V(T- 1 (y)) = vUh(\\y\\)-c,S- 1 



y 



WW) j ) 

= v[s~ x (||[|))-c + 7(lly||)=7(lly||) 

which finishes the proof of the first assertion. 
For n^4,5 and s > we define 

L(s) := sup \\DQ(x) 

V(x)=s 



and choose any class K, function a which is C°° and satisfies 

g 

a(s) < ——- for all s G (0, 1] 
L(s) 



Then the function h given by 

r 

h(r) = / a(s)ds 



o 



is smooth and of class /Coo. Note that this construction implies h(r) < ra(r) for all r > 0, 
hence h(r)/h'(r) < r. Thus 7 := /i" 1 is of class /Coo, smooth on (0, 00), and satisfies 

Differentiating T yields 

DT(x) = h'(V(x))Q(x) ■ DV(x) + /i(F(a:))DQ(a:). 

For rr — > the first term tends to since both h'(V(x)) = a(V(x)) and DV(x) tend to 0, 
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and the second tends to since for all x sufficiently close to the inequality 

h(V(x))\\DQ(x)\\ < a(V(x))\\DQ(x)\\ < ^-\\DQ(x)\\ < V(x) 



holds by construction of h. Thus DT(x) — > 0, as x — > 0, and consequently T G C 1 with 
DT(0) = 0, by a straightforward application of the mean value theorem, see e.g. [19, Chap. 
V, Theorem 3.2] and the fact that a function is continuously differentiable if all partial 
derivatives exist and are continuous. □ 



4 Main Results 



Using the coordinate transformation T we can now prove our main results. 

Theorem 2 Let n ^ 4,5 and consider any system (1) on W 1 which is UGAS (2). We 
suppose that the set D C W 1 is compact. Then, (1) can be transformed into a system (6) 
that is UGES (4). 

In particular, the constants in (4) can be chosen to be c = 1, A = 1. 



PROOF. Under our assumptions, by [20, Theorem 2.9, Remark 4.1]Q there exists a smooth 
function V : IT -> R for (1) such that 

L fd V(x) < - ai (\\x\\) (9) 



for some class /Coo function a%. Furthermore, there exist class /Coo functions c*2, «3 such that 

a 2 ( W) < V(x) < a 3 (\\x\\) . (10) 



Now let «4 be a C 1 function of class /Coo which is smooth on (0, oo) and satisfies o4(0) = 0, 
such that a 4 (a) < min{a, ct\ o (a)} for all a > 0. 

Such a function can be obtained e.g. by a slight modification of the construction in [26, Proof 
of Lemma 11]: Take a class /Coo function satisfying 6(a) < min{a, a\ o a^" 1 (a)} and which is 
smooth on (0, oo). Then 

aAa) = - / -dr 

71 J 1 + T 1 




* To be precise, the results in that reference make as a blanket assumption the hypothesis that / 
is locally Lipschitz, not merely continuous, at x = 0. However, as noted in e.g. [35], the Lipschitz 
condition at the origin is not used in the proofs. 
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has the desired properties. Thus we obtain 

L fd V(x) < -a 4 (V(x)). (11) 



Now define 



p{a) : = exp |— J 0:4 (t) 1 dr^ for a > 0, p(0) := 



Obviously p is smooth on (0, oo); furthermore p is of class /Coo and by [26, Lemma 12] p is 
a C 1 function on [0, oo) with p'(0) = 0. Thus defining 

W{x) :=p(V(x)) 



we obtain a C 1 Lyapunov function, which is smooth on M n \{0}, for which an easy calculation 
shows that 

exp ( — fy M aJr^dr) 



Applying Proposition 1 to W, using the class /Coo function 7 with 7(s)/7'(s) > s we obtain 
for each d E D and y 7^ 



Clearly the overall inequality also holds for y = so that we obtain 

%(t)f = 2(f(y(t),d(t)),y(t)) < -2\\y(t)f 



and hence ||y(t)|| 2 < e 2 *||y(0)|| 2 , i.e. the desired exponential estimate. □ 



Theorem 3 Let n 7^ 4,5 and suppose that the system (1) onW 1 is ISS (3) with some class 
/Cqo function a and some class ICC function j3. Then (1) is can be transformed into a system 
(6) that is ISES (5) with constants c = A = 1. 
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PROOF. By [32, Theorem there exists a C 1 function V which is smooth on IR n \ {0} 
and a class /Coo function x such that 

||s|| > X (||d||) =► L fd V(x)<- ai (\\x\\) 

for some class K,^ function a\. Furthermore, there exist class /Coo functions a 2 , a 3 such that 

a 2 (\\x\\)<V(x)<a 3 (\\x\\). 

As in the proof of Theorem 2 we find a function p which is class /Coo, C l , and smooth on 
R n \ {0}, such that W = p o V satisfies 

\\x\\ > X (\\d\\) => L fd W{x)<-W(x). 

Now Proposition 1 yields a parameter transformation T such that W(y) = PF(T _1 (?/)) = 
7(IMI) and t(s)/t / (s) > s. 

Now choose a class /Coo function 5 such that ||T _1 (?/)|| > 5(||y||) and define a = o x- 
Then a straightforward calculation yields 

||y||>fi(||d||) => ^W(y) < (12) 

Similar to the proof of Theorem 2 this implies 

||0(*,y,d(O)||<e-%|| 

as long as ||</>(t, ?/, > ( 5( su Po<r<< ll^( r )ll) which yields the desired estimate. □ 

Theorem 4 Consider the system (1) onl" being ISES (5) with some class /C^o function a 
and c = A = 1. T/ien there exists a homeomorphism R : R m — > IR m on the input space with 
-R(O) = 0, that is a diffeomorphism on lR m \ {0}, such that the the transformed system (7) 
satisfies the nonlinear Hoc estimate (8). 

PROOF. From (5) it is immediate that for any d(-) € T>, any x G M n , and any T > we 
have 

||x|| >e T a{ sup \\d(r)\\) => \\<f)(t,x,d(-))\\ < e-*\\x\\ for all te [0,T]. (13) 

0<t<T 

' As with the UGAS proof, it is easy to verify that the assumption that the right-hand side is 
Lipschitz at zero is never actually used in [32]. The possible non-uniqueness of trajectories does 
not affect the argument used in Lemma 2.12 in that paper, which reduces the problem to one of 
UGAS. 
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Now consider the function W(x) = \\x\\ 2 . Then (13) implies 

\\x\\ > e T a( sup \\d(r)\\) W(<j>(t,x,d(-))) < e~ 2t W(x) for all * G [0,T] . 

0<r<T 



In particular this estimate is valid for constant functions d(-) = d e D, thus the mean value 
theorem (observe W(<f>(0,x,d)) = W(x) = e~°W(x)) yields 

llxH > a(||d||) => L fd W(x) < -2W(x) < -W{x) . 

Now defining 

a(r) = sup (f(x,d),x) 

||x||<a(r),||d||<r 

we obtain a class K,^ function a with 

L fd W(x)<-W(x)+a(\\d\\). 



Without loss of generality (one could take a larger a) , we may assume a to be smooth on 
(0, oo), and thus 

R(d) := * m?d 



has the regularity properties as stated in the assertion. Now the transformation (7) yields 

LfW{x) < -W(x) + \\v\\ 2 . 

Integrating this equation along a trajectory x(-) gives 

i t 

W(x(t)) - W(x(0)) < -J W(x(s))ds + J \\v(s)\\ 2 ds 



which implies (8) since W(x) = \\x\\ 2 . □ 
5 Remarks 



Note that, in general, for our results to be true we cannot expect T to be diffeomorphic on 
the whole R™. Consider the simplest case where / does not depend on d and is differentiable 
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at the origin. If T were a diffeomorphism globally, then DT 1 (0) would be well-defined, 
which implies that 



DT{T-\y))f{T-\y)) = DT (0) D f (0) DT' 1 (0) 



y=0 



and so the linearizations in are similar; in particular, the dimension of center manifolds 
remains unchanged. 

Actually, if one wants the exponential decay to be e~*, even for linear systems one cannot 
obtain a diffeomorphism T. As an example, consider the one-dimensional system x = —x/2. 
Here one uses the change of variables y = T(x) given by T(x) = x 2 , x > 0, T(0) = and 
T(x) = —x 2 , x < to obtain y = —y. Note that T is C 1 with DT(0) = 0. The inverse of 
this T is given by T~ l (y) = ^/y, y > 0, T _1 (0) = and T~ 1 (y) = —y^y, y < which is 
smooth only away from the origin, though continuous globally. 

An example for the case of nontrivial center manifolds is given by the system x = —x 3 . Let us 
first note that for this system there is no transformation in the class we consider such that the 
transformed system is of the form y = —y. The reason for this is that we would have T(x) = 
y = —y = —T(x), so at least for x > V = T is a Lyapunov function with the property that 
V(x) = — V'(x)x 3 = — V(x). It is readily seen that the solutions of this differential equation 
(in x and V) are V c (x) = cexp^, for c G M. However, the image of [0, oo) under such 
V c yields a bounded set, so that these functions are no candidate for coordinate transforms 
on R. Nonetheless a coordinate transform according to our requirements can now be easily 
built: Take any /Cqo function a with a' > on (0, oo) so that with via the symmetrization 
a(—x) := a(x) we get a smooth function on R. Now define 

T(x) := a(x)Vi(x), x > 0, T(x) := — a(x)Vi(x), x < . 



Then for y ^ we have y = T(x) = — (1 + ~7~y )T(x) < —y, so that the transformed system 
decays at least exponentially with constants c = 1, A = 1. Again note that the requirement 
DT(0) = is vital, in fact all orders of derivatives vanish in 0. 

A basic ingredient of the proof of Theorem 2 is the construction of a Lyapunov function 
with the property V < —V. Actually, one may even, under restricted conditions, obtain 
the equality V — — V. It should be noted that already in [3] it is shown that for dynamical 
systems with globally asymptotically stable fixed point a continuous Lyapunov function with 
the property V((f>(t,x)) = e^V^x) exists, see also Chapter V.2 in [4]. Note, however that 
in these references only systems with trajectories defined on R are considered, which does 
not include the previous example. Indeed, if f(x, d) = f(x) is independent of d G D and the 
system x — f(x) is backward complete we can can also define a coordinate transformation 
based on a different W than the one used in the proof of Theorem 2: In this case the function 
W(x) = expt(x) with t(x) defined by V((f>(t(x),x)) = 1 is positive definite, proper, and 
satisfies L f W(x) = -W(x), thus W(<f>(t, x)) = W(x)-t. Since ^(l) = W'^l) we still find 
a diffeomorphism S as in the proof of Proposition 1. Deviating from this proof, instead of the 
gradient flow we now use the trajectories of the system, i.e. we define ir(x) = <f>(W(x) — 1, x) 
yielding W(n(x)) = W(x) — (W(x) — 1) = 1. Thus from ir we can construct T as in the 
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proof of Proposition 1, and obtain W(T -1 (y)) = \\y\\ 2 . Furthermore the definition of it 
implies that each trajectory {<p(t,x) \ t G M} is mapped onto the line {aS(7r(x)) \ a > 0} 
and consequently f(y) = —y, i.e. we obtain a transformation into the linear system y = —y. 
Note, however, that with this construction the coordinate transformation will in general only 
have the regularity of / (e.g. a homeomorphism if / is only C°), which is inevitable since 
it transforms / into a smooth map. Moreover, this construction cannot be generalized to 
systems with disturbances. 

Since we are not requiring that the inverse of a change of variables be itself a change of 
variables (because one may, and in fact does in our constructions, have DT(0) = 0, in which 
case T _1 is not differentiable at the origin), the way to define a notion of "equivalence" is by 
taking the transitive and symmetric closure of the relation given by such changes of variables. 
That is, we could say that system (1) is equivalent to a system (6) if there exist k G N and 
maps fo = f, fi, . . . , fk = f ■ IR n x D — > W 1 , all satisfying the assumptions on /, with the 
following properties: For each i — 0, . . . , k — 1 there exists a change of variables T as above 
such that fi(y, d) = DT{T- 1 (y))f m (T- 1 (y), d), where / = i, m = i + 1 or / = i + 1, m = i. 

Finally, regarding our notion of system transformation, note that even if /(0, d) ^ for some 
d G V for the original system (1), then under the assumption DT(0) = we have /(0, d) = 
for all d G T> for the transformed system. This implies that even if the original system had 
unique trajectories through zero, the transformed system cannot have this property. 

Acknowledgments: We thank David Angeli for suggestions regarding the remark on back- 
ward complete systems, as well as Uwe Helmke and Laurent Praly for many references to 
the literature. 
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